Dissipative hydrodynamics coupled to chiral fields 
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Using second-order dissipative hydrodynamics coupled self-consistently to the linear a model we 
study the 2 + 1 dimensional evolution of the fireball created in Au+Au relativistic collisions. We 
analyze the influence of the dynamics of the chiral fields on the charged-hadron elliptic flow V2 
and on the ratio Vi/{v2)^ for a temperature-independent as well as for a temperature-dependent 
viscosity-to-entropy ratio r]/ s calculated from the linearized Boltzmann equation in the relaxation 
time approximation. We find that «2 is not very sensitive to the coupling of chiral sources to the 
hydrodynamic evolution, but the temperature dependence of tj/s plays a much bigger role on this 
observable. On the other hand, the ratio V4,/{v2)'^ turns out to be much more sensitive than «2 to 
both the coupling of the chiral sources and the temperature dependence of r]/ s. 



I. INTRODUCTION 

One of the most surprising outcomes of the experi- 
ments conducted at the Relativistic Heavy Ion Collider 
(RHIC) is the discovery that high-energy collisions of 
heavy ions produce strongly interacting hadronic mat- 
ter (quark-gluon plasma - QGP) that evolves as a low- 
viscosity fluid [i|-|g|. This came as a surprise because 
early expectations, motivated mainly by the asymptotic 
freedom property of Quantum Chromodynamics (QCD), 
were that RHIC would create a gas-like system of weakly 
interacting quarks and gluons. 

The strongly interacting nature of the QGP is re- 
vealed through the description of measured momentum 
anisotropics via relativistic viscous hydrodynamics cal- 
culations. The momentum anisotropics are encoded 
in the Fourier moments V2, V3, V4, ■■■ of the mea- 
sured azimuthal distribution of particles, and are in- 
terpreted as being the translation to momentum space 
of the initial spatial eccentricity of non central colli- 
sions Q • A weakly-interacting, gas-like fluid would have 
no mechanism to induce such a translation into momen- 
tum anisotropics, but they can be produced if the parti- 
cles of the fluid are strongly interacting. 

In any attempt of describing experimental results of a 
heavy ion collision via a viscous hydrodynamics model, 
crucial physics input reflecting the properties of the flow- 
ing hadronic matter must be supplemented, like the equa- 
tion of state (EOS) and transport coefficients, as shear 
(77) and bulk {(,) viscosities. Although it is expected 
that for the strongly-coupled QGP 77 and ( will depend 
strongly on temperature - see below - usually in hydro- 
dynamic simulations temperature-independent values for 
these coefficients are assumed throughout the entire evo- 
lution. The impact of the temperature dependence of 
transport coefficients on momentum anisotropics that are 
obtained from hydrodynamic models has been recently 
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investigated in Refs. |8l-[l4l|. 

These transport coefflcients are, in principle, derivable 
from flrst-principles QCD calculations at finite temper- 
ature, but in practice, the calculations are extremely dif- 
ficult at strong coupling. There is a long history of com- 
putations of the EOS via large-scale numerical simula- 
tions of QCD on a lattice and modern calculations are 
achieving the accuracy needed for reliable use in heavy- 
ion physics - for a recent overview on the subject, see 
Ref. [15|. On the other hand, lattice QCD calculations 
of transport coefflcients are still in their infancy and only 
recently have results for shear and bulk viscosities been 
obtained [16l - ll8l |. Analytically, calculations have been 
performed within QCD for the shear and bulk viscosi- 
ties [1^ [20] at extremely high temperatures, where per- 
turbation theory in the coupling constant can be used 
on account of asymptotic freedom. Reliable calculations 
for these quantities can also be done at very low tem- 
peratures, where matter consists essentially of a very di- 
lute gas of pious. There is plenty of experimental infor- 
mation on elementary low-energy hadron-hadron inter- 
actions available that can be used to constrain detailed 
calculations that go beyond lowest-order in chiral per- 
turbation theory [2l|, [22]. On the other hand, in the 
region of intermediate temperatures, T ~ Aqcd, where 
one expects the deconfinement and chiral transitions to 
occur 23] , perturbative expansions (in the coupling con- 
stant or chiral) are not applicable. There is lack of experi- 
mental information on the appropriate degrees of freedom 
to use and phenomenological approaches are an alterna- 
tive, like the use of relativistic kinetic theory [24 30]. 

Assuming the validity of a relativistic viscous hydro- 
dynamical description, experimental extraction of trans- 
port coefficients should be possible in principle from com- 
parison with measured momentum anisotropy patterns. 
Moreover, as suggested in Ref. [3l[, such a comparison 
with experimental data would allow to pinpoint the lo- 
cation of the QCD phase transition or of a crossover from 
hadronic to quark-gluon matter. In this context, in ad- 
dition to the hydrodynamic degrees of freedom related 
to energy-momentum conservation, degrees of freedom 
associated with order parameters of broken continuous 



symmetries must be considered as well since they are 
all coupled to each other. Particularly important to the 
quest of determining possible signals of the chiral transi- 
tion is the coupling of chiral fields to the fiuid-dynamical 
modes. The present work is a first step towards the study 
of effects of such a coupling on momentum anisotropics 
within a viscous hydrodynamic description. 

Specifically, in the present paper we study, in the con- 
text of second-order dissipative relativistic hydrodynam- 
ics, the influence of the long-wavelength modes of chiral 
fields on the expansion of the fireball created in Au-|-Au 
collisions at y/s^N = 200 GeV. In particular, we aim 
at studying the effect of the coupled evolution of chiral 
degrees of freedom on the flow asymmetry characterized 
by V2 and V4^/{y2f', when finite viscosity is taken into 
account within a simple microscopic model for the chiral 
condensate. 

In our model the flowing matter consists of a fluid 
of quarks and antiquarks in local thermal equilibrium 
that evolves according to a second-order dissipative hy- 
drodynamic model [32h34| | with temperature-dependent 
speed of sound and transport coefficients. Starting from a 
high temperature state with approximately restored chi- 
ral symmetry, the system is evolved towards a state where 
the symmetry is spontaneously broken. The quark- 
antiquark fiuid interacts locally with the chiral fields of 
the linear a model (LSM) [33] ■ We compare the results 
obtained for hadronic observables when the chiral fields 
are included or not as sources in the hydrodynamic equa- 
tions. 

Similar models for the dynamics of quarks coupled to 
chiral fields were considered before. Mishustin and Scav- 
enius (36| constructed a relativistic mean field ideal fluid 
dynamical model based on the LSM, while Abada and 
Birse [37] used the Vlasov equation to describe the quarks 
coupled to the LSM. Both studies focused on the evo- 
lution of the chiral fields rather than on the influence 
of their dynamics on hadronic observables, which is the 
main focus here. Son [38| and Pujol and Davesne [3^ 
have developed modified hydrodynamic theories includ- 
ing additional (chiral) symmetry-breaking hydrodynamic 
variables, which were subsequently applied by Lallouet 
et al. [40| to the study of Bjorken flow in heavy-ion 
collisions. Later on, Paech et al. [4l], [i^ coupled 3-1-1 
ideal hydrodynamics to the LSM including fiuctuations 
of the chiral fields and studied the behavior of the az- 
imuthal momentum asymmetry at nonzero impact pa- 
rameter. More recently, Nahrgang and Bleicher [4^ in- 
cluded phenomenological dissipation and noise terms in 
the equation of motion of the a field of the LSM coupled 
to ideal hydrodynamics. They focused on sigma fluctu- 
ations and found that at the flrst order phase transition 
they lead to an increase in the intensity of sigma excita- 
tions. In a subsequent work, Nahrgang, Leupold, Herold 
and Bleicher ^J] employed the two-particle irreducible 
(2PI) effective action formalism to set up an approxi- 
mation scheme to derive explicit formulas for the dissi- 
pation kernel and the noise correlation function within 



the LSM. While the approximations employed deserve 
careful scrutiny, the formulas derived within the men- 
tioned approximation scheme show interesting features 
of dissipation and noise on the relaxation dynamics of 
the chiral dynamics and point to the importance of in- 
cluding reheating of the quark fiuid by energy dissipation 
from the chiral fields, as shown in the very recent papers 
by Nahrgang, Leupold and Bleicher Ii45il and Nahrgang, 
Bleicher, Leupold, and Mishustin [4a|. The work of 
Plumari et al. [43] is more in the spirit of ours in the 
present paper, in that they have investigated the role 
of the chiral transition on V2- Specifically, the authors 
have solved numerically the Boltzmann- Vlasov transport 
equations including both two-body collisions and the chi- 
ral phase transition as given by the Nambu-Jona-Lasinio 
model. For conditions prevailing in Au -I- Au collisions at 
y/sNj\j = 200 GeV, the authors find a sizable suppression 
of V2 due to the attractive nature of the chiral fleld dy- 
namics. They also found that, if rj/s is kept fixed, V2 does 
not depend on the details of the collisional and/or field 
dynamics and in particular it is not affected significantly 
by the chiral phase transition. 



At this point it is appropriate to state the limited scope 
of our work. First of all, while we use a viscous hydro 
model coupled to a workable chiral model, we do not 
intend to extract values of transport coefficients by com- 
paring our results to V2 data. Actually, extraction of 
transport coefficients by matching the output of hydro- 
dynamic and/or kinetic simulations to data is not at all 
straightforward, as discussed e.g. in Refs. [ol. 1481 - 1501 ] . In 
particular, to extract rj/s from V2 measurements, a pos- 
sible route is to first constrain initial Ti and final Tj 
temperatures (that are used to run the coupled set of 
hydro-chiral equations) by matching calculated spectra 
to data and, afterwards, extract r//s from V2 data - see 
e.g. Refs. [43,[51|. In addition, although we use a semi- 
realistic initial condition (we use a Glauber smooth initial 
condition), a fluctuating, lumpy initial condition leads to 
lower elliptic flow along with other interesting effects |52| - 
1541 ] that certainly have an impact on the extraction of 
transport coefficients. 



The paper is organized as follows. In Section |TI] we 
describe the coupling of the chiral fields to the 2-1-1 
hydrodynamic model. In Section Hill we describe the ini- 
tial conditions and the freeze-out prescription used in 
the simulations. We then show and discuss our results 
in Section IIVI and conclude in Section |Vl In Appendix 
\X\ we briefly overview the LSM and the calculation of 
the speed of sound and of the shear viscosity from the 
linearized Boltzmann equation in the relaxation time ap- 
proximation. In Appendix [B] we discuss the dependence 
of our results on the cut-off that we must impose on the 
value of rj/s due to the breakdown of viscous hydrody- 
namics when nonequilibrium effects start to dominate. 



II. COUPLED CHIRAL-HYDRO DYNAMICS 

In this section we describe the model of couphng chi- 
ral fields to hydrodynamical variables. The evolution of 
the chiral fields is described by the LSM [35| and hy- 



drodynamics is described by dissipative hydrodynamic 
equations |32h34| |. This hydrodynamic formalism goes 
beyond the well-known Israel-Stewart theory [55l - l57l | in 
that it includes all second-order terms in velocity gra- 
dients that can appear in the stress-energy tensor of a 
conformal fluid - see e.g. Ref. J58J. This hydro model 
has been applied to model the expansion of the QGP in 
Refs. m and [Hj. 

In order to set notation and make the paper self- 
contained, we start with the LSM model. The La- 
grangian density of the model, given in terms of quark 
q = (u, d) and chiral (j)a = (cr, n) fields, is written as (we 
use the signature for the metric g^i^ = (-I-, — , — , — ) ) 



with 



A-^ *^q ~r A-'rf), 



^q^q [il'^df, ~ g{(T + j5T- 7?)] q, 



(1) 



(2) 



p = p{(j), T) are the energy density and pressure (by def- 
inition) at equilibrium at local temperature T ~ their 
explicit expressions are given in Appendix |Al The con- 
tribution from the chiral fields to the energy-momentum 
tensor is given by 
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a>,-g^^(/:^). 



(7) 



where (• ■ ■ ) means local thermal average. 

The hydrodynamic equations of motion are obtained 
from conservation of the energy-momentum tensor, 
D^T'^'^ — 0, where D^ is the geometric covariant deriva- 
tive. We employ Milne coordinates (appropriate for a 
(2 + 1) fiow dynamics) defined by proper time r = 
y/t"^ — z^ and rapidity ^ — arctanh(2;/t). The hydro- 
dynamic velocity is u = (m'^,m^,u^,0) and the metric 
tensor reads g^^ — (—1, 1, Ij^^). In the transverse plane 
we use Cartesian coordinates x^ — {x,y), so that the 
only non-vanishing Christoffel symbols are ri,^, = r and 

r^^ — 1/t. Since we assume boost-invariance all quanti- 
ties are independent of rapidity. The conservation equa- 
tions then read 



and 



where 



^0=2 (d^'^d'^'y + d^^ ■ d^^) - U{<J, tt). 



(3) 



C/(a,^) 



U{<f>) 
A2 , 



- hqa - Uq 



'Y-hqa-Uo, (4) 



is the potential which exhibits chiral symmetry break- 
ing - with 0^ = cr^ -I- TT^ . At the level of approximation 
we work in the present paper, parameters are fixed as 
follows. The vacuum expectation values of the conden- 
sates are taken (a) — /^r and (tt) — 0, with the pion 
decay constant /^r — 93 MeV. Also, the partially con- 
served axial-vector current relation yields hq — /tt^z^ 



m: 



\/y 



with m^ = 138 MeV. This leads to u^ = fl 
A mass nia — \/2\^f'^ -\- m^ ~ 600 MeV is obtained if 
A^ = 20. The constant Uq is chosen such that the poten- 
tial energy vanishes in the ground state. 

Following Paech et al. [4l|, |42|, we split the energy- 
momentum tensor as 



rpf_il' rpf_iU 



T 



ifiiy 



(5) 



As mentioned previously, the quarks and antiquarks are 
assumed to constitute a heat bath in local thermal equi- 
librium and their dynamical evolution will be determined 
by viscous relativistic hydrodynamics, so that: 



q 



{e+p)u''u-' -pg"" +n'"', 



(6) 



where u^ is the four- velocity normalized as u^u'^ = 1, 
n^^" is the viscous shear-tensor, and e = e(0, T) and 



+ g {psDa + pps ■ Dtt) , (8) 

(e + p) Du' = cl (g'^djE - u'u'^do^e) - IV^DpIi"''^ 

-f-5(p,VV-f pp, • W) , (9) 

where A^" — g^^'^ — u^^u'^ is a projector (orthogonal to the 
fluid velocity), Cg is the speed of sound, and ps = (qq) 
and Pps — (qj^Tq) are the local thermal averages of the 
scalar and pseudo-vector chiral densities - their explicit 
expressions are given in Appendix [Xj It is seen that the 
chiral densities act as sources for the evolution of the 
hydrodynamic variables e and u'^. Here, D = UpD'^ is 
the comoving time derivative, V^ = A^qZ?" is the spatial 
gradient, a^i, is the shear tensor: 



a'"' = V^'^u'^^ 



(10) 



The angular braces around Lorentz indices denotes the 
spatial, symmetric and traceless part of a tensor, i.e. if 
A^" is a tensor then A^^^'^'^ means: 

The evolution of the shear tensor depends only indirectly 
on the chiral fields and it is given by 
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(12) 



where i] is the shear viscosity and (r,r,Ai) are second- 
order transport coefficients. We note that Israel-Stewart 
formahsm is recovered from these equations when Ai = 0. 
In the Milne coordinates, the classical equations of mo- 
tion for the chiral fields evolving in the background of the 
quark fluid read 



T J 5a 



di + -dr - df I vr 



SU_ 
Sir 



= -gps , 



= -gpps 



(13) 



(14) 



where ps and ppg are the scalar and pseudo-scalar chiral 
condensates defined above. 

In this way, the evolution of the chiral fields affects 
the evolution of the quark fluid through the sources in 
the energy-momentum conservation equations; in turn, 
the quark fluid affects the evolution of the chiral fields 
through the densities ps and pps- 



III. SOLVING THE EQUATIONS 

We solve the hydrodynamic equations given in Eqs. (|5]) 
and dH) using the publicly available code of Luzum 
and Romatschke [i^. We need to supply the locally 
temperature-dependent transport coefficients r], r^ and 
Ai, and speed of sound Cg- In addition, we need to spec- 
ify initial conditions for the energy density e, velocities ui 
and U2, and the components of the shear tensor. To solve 
the equations of motion for the chiral fields, Eqs. P^ and 
(HH) , we use a simple finite-difference scheme both in the 
proper-time and space variables. Here also, one needs to 
provide initial conditions for the fields and their deriva- 
tives. 

In all calculations performed in the present paper, we 
use a 13 fm X 13 fm transverse plane. In the next subsec- 
tions we discuss the values used for the input constants 
and initial conditions. 



A. Transport coefficients and the speed of sound 

As a prototypical example of a relativistic fluid, we 
will use r^ = 2(2 - ln2)?7/sT and Ai = ?y/27rT, that 
correspond - in a gradient expansion at second-order - 
to a supersymmetric Super- Yang-Mills plasma |3a - l34 |. 
However, as mentioned in the Introduction, we will not 
use c^ = 1/3 and rj/s — 1/4tt appropriate to the Super- 
Yang-Mills plasma. The temperature dependence of the 
speed of sound is obtained within the LSM and that of 
J] is obtained using the linearized Boltzmann equation 
in the relaxation time approximation within the chiral 
model of Ref. [2^ . We discuss these in the Appendix |^ 

We note that the value that we use for the relaxation 
time Ttt = 2(2 — ln2)77/sT is not the exact value of the 
Super- Yang-Mills plasma [5^, |60| . Rather, it is the 



value obtained through a second-order gradient expan- 
sion, that is equivalent to a Taylor expansion of the re- 
tarded Green's function that linearly relates a dissipative 
current and a thermodynamical force. The expression 
for r^ derived in this way is in general different from the 
one obtained directly from the first pole of the retarded 
Green's function. Moreover, as shown in Refs. |59l. [60|. 
the true dynamics of the long-wavelength and low fre- 
quency modes of strongly coupled theories is not gener- 
ally described by relaxation-type equations, i.e., the evo- 
lution equation for the shear tensor includes terms with 
second time derivatives of the thermodynamic forces - 
see also Ref. 6l|. 

In this connection, it is important to remark that it 
was shown in previous hydrodynamic simulations (with- 
out chiral fields) [4^, [51j that, for values of rj/s < 0.2, 
the results for the hadronic observables do not depend 
strongly on the choices for the second-order transport 
coefficients, provided that Ai 7^ 0. The infiuence of terms 
with second time derivatives of the thermodynamic forces 
in the evolution of the shear tensor of the QGP is still 
an open issue that surely deserves further investigation. 
In this work we will neglect such terms and use instead 
the relaxation-type equation given in Eq. P^ . that is 
obtained from a second-order gradient expansion. 



B. Initial conditions 

For the initial transverse velocity and shear tensor we 
use M^ = w** = 0, which implies vanishing initial vor- 
ticity, and (n^^,n^2', H^'^') = 0. The initialization time 
is set to To = 1 fm/c. It has been shown before that 
the evolution of the shear tensor H'^'' is quite insensitive 
to the initialization values, the difference being appre- 
ciable only at very early times (see e.g. Ref. [50|). We 
have verified that the elliptic flow shows very little sen- 
sitivity to the initialization of the shear tensor as well. 
The initial temperature at the center of the flreball is 
set to Ti — 333 MeV. This value for Ti has been used 
in previous hydrodynamic simulations and leads to kaon 
multiplicity and (pt) that are in good agreement with 
RHIC data ^^. 

The initial energy density proflle is calculated using 
Glauber's model, in which for a given impact parameter 
b we have 

eiTo,x,y,b) = CaTAix + b/2,y)TA{x-b/2,y) , (15) 

where C is a constant chosen such that e(ro, 0,0,0) cor- 
responds to a given initialization temperature Tq via the 
EOS, the cross-section a is taken to be cr = 40 mb, and 
Ta is the nuclear thickness function given by 



TA{x,y)=l 5A{x,y,z)dz, (16) 

J —00 

with 5a is the density distribution for the gold nucleus, 



taken to be of a Woods-Saxon form 

So 



6Aix,y,z) 



l + exp[(|f|-i?o)/x] ' 



(17) 



where x — {x, y, z), Rq = 6.4 fm and x = 0.54 fm. The 
parameter Jq is chosen such that / d^x 5yi(x) = 197, as 
appropriate for Au nuclei. 

As a reasonable ansatz for the initial condition of the 
chiral fields we use 



Tf{To,r) = and cr(ro,r) = /^ 



1 



-('■/'■o)' 



(18) 



with r^ 



y^ and rg = 9 fm. In this way, the 



chiral condensate nearly vanishes at the center where the 
temperature of the fluid is larger and interpolates to f.^ 
where the temperature is lower. 



C. Freeze-out scheme 

In our simulations we use the isothermal Cooper-Frye 
freeze-out prescription [62| in which the conversion from 
hydrodynamic to particle degrees of freedom takes place 
in a three-dimensional hypersurface. The spectrum for a 
single on-shell particle with momentum p^ = (E^p) and 
degeneracy d is 



E 



dN 



d^p (27r)3 



p^ dE^ f{x^,pn 



(19) 



where dE'^ is the normal vector on the hypersurface. 
The non-equilibrium distribution function / is given by 
Grad's quadratic ansatz |63l |: 



f{x^,pn = 



/o(a:^P^) 

-t- fo{x^,pn [1 T /o(x^p^)] /ff/7e) 



1 



PmP-H''" 



2T^ip + e) 



exp 



-Pp' 



T 



(20) 



where /o is the Fermi-Dirac distribution. The approxi- 
mation in the third line holds when p >> T, and it is 
used in our simulations. The systematic error of this ap- 
proximation is very small at low tranverse momentum 
Pt < 2.5 GeV, so we do not expect our results to have 
a significant error coming from this approximation (see 
Ref. J&]). 

We calculate the spectra for particle resonances with 
masses up to 2 GeV and then determine the spectra of 
stable particles including feed-down contributions. For 
this last step we use the AZHYDRO package [6J-[6g. 
In the present paper we will focus on the elliptic flow 
coefficient V2 and also on Vi at central rapidity, which 
are related to the particle spectra (including feed-down 
contributions) by 



E 



.dN 
d^p 



vq 



{pT,b)[l + J2^MPT,b)cos{nip)] , (21) 



with (f = a,Tcta,n{py/px) and pr = {p^ + PyY^^- For 
the freeze-out temperature we take Tf = 130 MeV. This 
value for Tf is slightly smaller than the one usually em- 
ployed in viscous hydrodynamic simulations of Au-|-Au 
collisions of Tf ^ 140 MeV, but allows us to study a 
broader range of temperatures near our value of Tc in 
order to determine the influence of chiral fields dynamics 
on hadronic observables. 



IV. RESULTS 

We now go over to discuss our results. We start with a 
very important input to the hydro equations, the square 
of the speed of sound c^ = dp/de. In Fig. [T] we show 
cl as predicted by the LSM - its explicit expression is 
given in Appendix |X] - for three values of the coupling 
constant, namely g = 3.0, 3.2, 3.4, which correspond to 
a smooth crossover. The reason to consider these values 
for g is that recent works have shown that the crossover 
phase transition without a thermodynamic region where 
the sound velocity drops to zero leads to a faster time 
development of the system and helps to reproduce RHIG 
data with hydrodynamic simulations |67h69| . Moreover, 
Lattice QCD calculations also favor a smooth crossover 
with a relatively soft dip in c? over a first-order phase 
transition (see e.g. Refs. |7C)I - [74| '). In our case, the first- 
order transition is obtained when g = 3.8. As it can 
be seen from Fig. [2 lowering the value of g leads to 
a softening of the crossover and a raise in the critical 
temperature Tc- For large temperatures the conformal 
limit c^ = 1/3 is reached, while for low values c^ goes to 
zero. For g ^ 3.2 the behavior of c^ with temperature 
is in qualitative agreement with that obtained in Lattice 
QCD calculations [7Cll - l74| and with that favored by hy- 
drodynamic simulations based on different interpolation 
schemes used to join Lattice QCD and hadron resonance 
gas equations of state (see e.g. Refs. [67H6g |). For g > 3.4 
or g < 3, the dip in c^ at Tc becomes very sharp or fades 
away, respectively, resulting in a temperature dependence 
that is in disagreement with the one obtained in Lattice 
QCD calculations. 

The critical temperature for these values of the LSM 
parameters is Tc ~ 150 MeV, which is smaller than the 
value ~ 170 MeV obtained from Lattice QCD simulations 
[70l - [74| . Within the LSM it is not possible to obtain much 
larger values for Tc without going to very low values of 
g. These lower values of g correspond to a very soft 
crossover with an almost inappreciable drop in cj. that is 
not in agreement with the results the Lattice QCD. The 
relatively low value for Tc used here does not constitute 
a serious limitation to our aim of obtaining a qualitative 
understanding of the impact of the evolution of the chiral 
fields on hadronic observables, since we expect the trends 
obtained in this paper to hold in more realistic scenarios. 



Another input to the hydrodynamic simulation is the 
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FIG. 1: (Color online) Square of the sound speed c^ of the 
LSM as a function of temperature, for three values of the 
chiral coupling: g = 3, 3.2, 3.4. 



ratio r]/s as a function of temperature, which we calcu- 
late in the LSM from the linearized Boltzmann equation 
in the relaxation time approximation (details are given 
in Appendix El) . Fig. [5] shows ij/s as a function of tem- 
perature in a region close to Tc, for g — 3,3.2,3.4. It is 
seen that rj/s decreases significantly when approaching 
Tc from below, and remains almost constant and rather 
small at T > Tc. This behavior is qualitatively similar to 
that found in other approaches such as the Boltzmann- 
hydrodynamics hybrid approach [73, the NJL model [2M| 
and Chiral Perturbation Theory [76|, although, as al- 
ready indicated, the value of Tc obtained here is smaller. 
It is also seen that when the value of g increases the vari- 
ation of r]/s with temperature is more abrupt, decreasing 
faster at T < T,, and thus remaining small for a larger 
range of temperatures. 
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FIG. 2: (Color online) rj/s of the LSM as a function of temper- 
ature, for three values of the chiral coupling: g — 3, 3.2, 3.4. 

At this stage it is appropriate to note an important 
point concerning the rise of rj/s with decreasing tem- 



perature in the low temperature stage. By construc- 
tion, any hydrodynamic description of matter necessar- 
ily breaks down when the viscosity is large enough that 
nonequilibrium effects are as important, or even more 
important than the equilibrium ones. For example in 
the context of Grad's quadratic ansatz that is usually 
employed at freeze-out, this occurs when the nonequilib- 
rium correction to the distribution function is compara- 
ble - or exceeds - the equilibrium one. In the context of 
relativistic heavy-ion collisions numerous hydrodynamic 
simulations (usually with temperature-independent rj/s) 
have shown that large enough means rj/s ~ 0.3 — 0.5 
O [H, US III [zS-il. For larger values of rj/s it is 
not possible to match the predictions of hydrodynamic 
models to RHIC data. Moreover, it has been shown by 
direct comparison to results obtained with Boltzmann 
equation that for matter created at RHIC, viscous hy- 
drodynamics fails to reproduce the results of Boltzmann 
equation if rj/s > 0.2 |8Ij . Beyond this value, an appro- 
priate treatment of the dynamics of this matter should 
switch smoothly from a viscous hydrodynamic descrip- 
tion to a kinetic one, in which the issue of breakdown 
due to large viscosity (i.e. dissipation) does not arise. 
For recent developments and applications of this hybrid 
description see Refs. [l0ll75|. 

Since the model used here is purely hydrodynamic, a 
cut-off for the value of rj/s must be imposed in order to 
avoid the breakdown of the fluid description. The value 
of rj/s at which it is sensible to impose this cut-off is 
constrained by the comparison of results obtained from 
hydrodynamic and kinetic simulations to data. For this 
reason, in our simulations we set rj/s < 0.4, which corre- 
sponds to T ^ 440 MeV for the three values of the chiral 
coupling, namely g = 3, 3.2, 3.4, the values we consider 
in the present paper. We emphasize that in all of our 
simulations with g = 3,3.2,3.4 the average value of rj/s 
throughout the hydrodynamic evolution is ~ 0.11, and 
thus it remains well below the value at which viscous 
hydrodynamics breaks down {rj/s ^ 0.3 — 0.5). 

In Appendix [B] we analyse the dependence of our re- 
sults on the choice of different cut-off values for tj/s. 
We find that, although the values of V2 and V4/{v2)'^ do 
change, the differences between the cases including or not 
the chiral fields as sources remain essentially the same, 
and thus our main conclusions are not sensitive to the 
value imposed for the cut-off. 

Knowing c^ and rj/s as functions of the tempera- 
ture, we can now go over to discuss the momentum 
anisotropics obtained from the simulations. Fig. |3] shows 
the charged-hadron elliptic flow V2 calculated with ei- 
ther the temperature-dependent rj/s of the LSM or a 
temperature-independent rj/s = 0.11, which corresponds 
to the average of rj/s throughout the hydrodynamic evo- 
lution over the temperature interval from T^ to Tf. In 
both cases, we set g = 3.2 and compare the results ob- 
tained by taking or not taking into account the source 
terms in the hydrodynamic equations. 

It is seen from Fig. [3] that both for a temperature- 




FIG. 3: (Color online) Elliptic flow of charged-hadrons calcu- 
lated taking or not taking into account the source terms in the 
hydrodynamic equations, for either a temperature-dependent 
or a temperature-independent -q/s which is the averaged value 
throughout the evolution of the fireball. The value of the cou- 
pling constant is g = 3.2. 



dependent and a temperature-independent 77/s, the el- 
liptic flovif does not depend strongly on the chiral sources. 
For a temperature-independent 77/s, V2 turns out to 
be practically independent of the dynamics of the chi- 
ral sources for pT < ?> GeV. In contrast, for the 
temperature-dependent 77/s there are small, but visible 
differences in the elliptic flow calculated with or without 
the chiral sources. Specifically, V2 is slightly larger and 
reaches a maximum at pt ^ 2 GeV when the chiral fields 
are taken into account as sources for the hydrodynamic 
variables. 

Comparing the behavior of V2 for temperature- 
independent or temperature-dependent rj/s, it is seen 
that a temperature-dependent 77/s leads to a smaller 
elliptic flow. This fact has already been seen in the 
studies of Refs. [g, [3, \lM by performing hydrodynamic 
simulations (not coupled to chiral fields) with different 
parameterizations for temperature-dependent 77/s (usu- 
ally a step-function dependence or some slight variation), 
and in Ref. [75| from a Boltzmann-hydrodynamics hybrid 
approach. See also Refs. [10, ll^ for related studies. In 
particular, it was found in Ref. Q that V2 at Au-|-Au col- 
lisions at the largest RHIC energy is dominated by the 
shear viscosity in the hadronic phase and largely insen- 
sitive to the viscosity of the QGP. A similar conclusion 
was reached in Refs. [^, [il|. It is interesting to note 
that, according to recent viscous hydrodynamic simula- 
tions with temperature dependent 7//s, at LHC energies 
the situation is exactly the opposite, i.e. the elliptic flow 
becomes sensitive to the QGP viscosity and insensitive 
to the hadronic viscosity [a, Q • As noted in Ref. Q , the 
fact V2 depends rather weakly on the temperature depen- 
dence of ?7/s poses serious challenges to the precise ex- 
traction of this ratio from this hadronic observable. For 
this reason, and also due to the uncertainty in the ini- 



tial conditions in heavy-ion collisions at RHIC and LHC, 
it has now become important to calculate higher-order 
Fourier moments V3, ■ ■ ■ in viscous hydrodynamic simu- 
lations, which may provide further crucial constraints on 
77/s and on models used to calculate QGP initial condi- 
tions - see e.g. Refs. [2i,|8i-[83. 

Returning to our results in Fig. |31 it seems clear the 
difficulties in extracting (an average value of) ?7/s from 
data on V2. The results show that uncertainties associ- 
ated with the dependence of 77/s on temperature lead to 
appreciable changes in the curve of V2 versus pt- Such 
an uncertainty should be added to the theoretical uncer- 
tainty that comes e.g. from the initial conditions (for 
example using Color Glass Condensate or Glauber initial 
conditions), the freeze-out process (as issues concerning 
Grad's quadratic ansatz), as well as from other sources, 
that according to recent studies add up to an overall un- 
certainty which can be roughly estimated in 0.1 |48l.l5l|. 

As mentioned in the Introduction, the fact that V2 de- 
pends somewhat on the evolution of the chiral fields in 
the case of a temperature-dependent 77/s has already 
been seen in the study by Plumari et al [47| within 
Boltzmann-Vlasov simulations on the NJL model. Our 
results show that this also happens in a purely hydrody- 
namic simulation as well. 
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FIG. 4: (Color online) Elliptic flow of charged-hadrons cal- 
culated taking or not taking into account the source terms 
in the hydrodynamic equations, for g = 3, 3.2, 3.4 and a 
temperature-dependent rj/s. 



Since our approach is phenomenological, it is impor- 
tant to determine the sensitivity of our conclusions re- 
garding the infiuence of the chiral fields on V2 to the value 
of the coupling constant g. Fig. 2] shows the elliptic flow 
of charged-hadrons calculated taking or not taking into 
account the source terms in the hydrodynamic equations, 
for g = 3, 3.2, 3.4 and a temperature-dependent 7;/s. It 
is seen that when the value of g increases, V2 decreases. 
This is due to the fact that for larger values of g, c^ near 
the transition region becomes smaller which means that 
pressure gradients are converted into flow less efficiently, 
thus resulting in lower values of V2 ■ The difference in the 
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values of V2 calculated including or not the chiral fields 
as sources is seen to become smaller with increasing g, 
which supports the conclusion reached earlier that the 
influence of the chiral fields on V2 is rather small. 
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FIG. 5: (Color online) V4/{v2)'^ for charged-hadrons cal- 
culated taking or not taking into account the chiral fields 
as sources for the hydrodynamic equations, for either a 
temperature- dependent or a temperature-independent rj/s. 
The value of the coupling constant is g = 3.2. 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 

0.05 



-0.05 



•' >«jS'"""""., 




g=3.4 w/sources 

g=3.4 no sources * 

=3.2 w/sources 

g=3.2 no sources « 

g=3 w/sources 

=3 no sources 



0.5 



1.5 2 

Pt [GeVl 



2.5 



3.5 



FIG. 6: (Color online) 1)4/(1)2)^ for charged-hadrons calcu- 
lated taking or not taking into account the chiral fields as 
sources for the hydrodynamic equations, for g = 3, 3.2, 3.4 
and a temperature-dependent 77/s. 

Another interesting observable to analyse is V4/{v2)'^, 
which has only recently been theoretically investigated. 
See in particular Refs. [83, |80| for a detailed account of 
the physics (and caveats - see below) involved in this 
observable within dissipative hydrodynamics. See also 
Ref. \47l for a calculation of Vi/{v2)^ within Boltzmann- 
Vlasov formalism. Fig. [5] shows W4/(w2)^ for the same 
cases as in Fig. [3l In contrast to what happens with 
^2, it is seen that f4/(f2)^ is significantly affected by 
the dynamics of the chiral fields even at low transverse 



momentum and for a temperature-independent 77/s. In 
particular, the values of f4/(u2)^ are larger when the chi- 
ral fields are included as sources in the hydrodynamic 
equations, and the behavior with px is somewhat differ- 
ent too. It is also seen that a temperature-independent 
T]/s leads to larger V4/{v2)^- 

At this point, a comment on the possible extraction of 
rj/s from f4/(f2)^ data by matching to hydrodynamic 
simulations is in order. It is known that the depen- 
dence of ^4/(^2)^ on tranvserse momentum that is ob- 
tained from viscous hydrodynamics does not agree with 
the almost constant value of V4/{v2)'^ ~ 1 measured at 
RHIC (see Refs. [83, [S^]). Ideal hydrodynamics predicts 
'^'4/(^2)^ ~ 0.5, but quite surprisingly viscous hydro- 
dynamics yields a strongly dependent ratio, as shown, 
for example, in Fig. [S] A possible explanation for this 
discrepancy has been put forward by Luzum and OUi- 
trault 88', 89], who realized that Grad's quadratic ansatz 
for the nonequilibrium correction 6f to the particle distri- 
bution fuction may not be valid for the freeze-out process 
in heavy-ion collisions. See also Refs. |90l - l93j for recent 
work on the reliability of Grad's ansatz. By performing 
viscous hydrodynamic simulations with different depen- 
dencies of 5f on Pt, the authors of Refs. [8^, [83] have 
found that RHIC data on W4/(w2)^ favors a momentum 
dependence between linear and quadratic. In spite of 
this, we have chosen to show results on W4/(f2)^ to em- 
phasize that, although the values for this observable are 
very different from those measured at RHIC (because we 
are using Grad's quadratic ansatz at freeze-out), they are 
much more sensitive than V2 to the evolution of the chiral 
fields. We expect this feature of ^4/(^2)^ to hold when 
other forms for Sf are used in the simulations, although 
further work is needed to confirm this expectation. 

As with V2, it is important to determine the depen- 
dence of our results and conclusions on the value of g. 
Fig. [5] shows W4/(w2)^ calculated including or not the 
sources in the hydrodynamic equations, for g = 3, 3.2, 3.4 
and a temperature-dependent rj/s. It is seen that the 
difference in the values of W4/(w2)^ calculated including 
or not the chiral fields as sources becomes smaller with 
increasing g, but is still relatively large at g — 3.4. More- 
over, it is seen that the behavior of this ratio with px does 
not depend strongly on the value of g. These results ob- 
tained with different values of g support the conclusion 
mentioned before that V4/(v2)'^ is more sensitive than V2 
to the influence of the chiral fields on the evolution of the 
quark fluid. 



V. CONCLUSIONS 

We have studied the evolution of the flreball created at 
RHIC within second-order viscous hydrodynamics cou- 
pled self-consistently to the LSM, focusing on the impact 
of the dynamics of the chiral fields on the hadronic ob- 
servables V2 and V4/{v2)'^. We have compared the results 
obtained when the chiral fields are included or not as 



sources in the hydrodynamic equations. The compar- 
isons were made using both temperature-independent 
and temperature-dependent 77/s. The temperature de- 
pendence of 77/s was calculated in the LSM using the 
linearized Boltzmann equation. 

We have found that the values of V2 do not depend 
strongly on the evolution of the chiral fields. Specifically, 
for a temperature-independent 77/s this dependence is 
negligible for px < 3 GeV, while for a temperature- 
dependent r]/s it is appreciable but still small even at 
small pt- We have also found that the ratio Vi/{v2)^ is 
much more sensitive to the dynamics of the chiral fields, 
being larger when these fields are taken into account as 
sources in the hydrodynamic equations, in both situa- 
tions of the temperature dependence of 77/s. 

In line with the results of Refs. |8l-[l(]||. our results show 
that not knowing precisely the temperature-dependence 
of 77/s leads to further uncertainties in attempts of ex- 
tracting this ratio from data on V2, in addition to the 
uncertainties that stem from the initial conditions and 
the freeze-out process, among others sources. It is worth 
noting that despite the coupling of chiral sources to the 
hydrodynamic evolution would add further uncertainties, 
they are not very big. 

The model used in this work leaves room for improve- 
ments in different directions. Probably, the most impor- 
tant ones for our analysis are including bulk viscosity in 
the hydrodynamic equations and fluctuations of the chi- 
ral fields. This latter effect would act as noise sources 
in the classical equations of motion. We believe that the 
qualitative trends and the general conclusions extracted 
from our results will hold when these effects are taken 
into account. Work is in progress where these two as- 
pects are taken into account in a viscous hydrodynamic 
simulation and their impact on observables will be re- 
ported in a forthcoming publication. 



Appendix A: Linear a model 

In this Appendix we will briefly review the relevant 
aspects of linear a model for the present paper, as well 
as the calculation of 77 from the linearized Boltzmann 
equation in the relaxation time approximation. 

As mentioned before, as an effective theory of the chi- 
ral symmetry breaking dynamics we consider the linear 
a model coupled to two flavors of constitutive quarks. 
The Lagrangian density of the coupled system is given 
in Eq. ([T]). The quark fluid is considered as a ther- 
mal bath for the chiral fleld, and therefore it can be 
integrated out to obtain the effective potential Ve for 
the chiral flelds in the presence of that bath of quarks. 
In order to calculate the equilibrium pressure density 
75(0, T) = —Ve(<t>, T) -!-[/((/>), we use the one-loop effective 
potential [4ll.l4a|: 



v;(0,r)-c/(0) 



dqT 



_fp_ 

(2^ 



In 



{l + e'^/^) , (Al) 



24 is the color-spin-isospin-baryon charge 
if the quarks and the energy E = {p^+m'^Y^'^ 
^0^. The equilibrium energy density is just 



e(0,r) 



d^p 

(2^ 



hE 



(A2) 



From p and e, the square of the speed of sound and the 
entropy density follow: c^ = dp/de and s = dp/dT . 

The local thermal averages of the scalar and pseudo- 
vector chiral densities ps = {qq) and pps = {q^y^rq) enter 
as sources in the hydrodynamic equations. Explicitly, 
they are given by: 



Ps 
Pps 



gadg 



d^p 1 

(27r)3 E 



/o, 



-^ I d'p 1 ^ 



(A3) 
(A4) 



The LSM exhibits a flrst-order phase transition, a 
crossover and a critical end point, depending on the 
value of the chiral coupling constant g. The crossover 
phase transition is not a genuine phase transition since 
all thermodynamic functions change smoothly with tem- 
perature. However, such changes may be quite sudden in 
a narrow temperature range. As indicated in Section HVl 
throughout this work we set g — 3, g = 3.2 and g = 3.4, 
which correspond to a smooth crossover and lead to a 
temperature-dependent sound speed that resembles that 
of Lattice QCD. 



^ 0.6 




FIG. 7: (Color online) (e - Sp)/T'^ of the LSM as a function 
of temperature for g — 3, 3.2, 3.4. 

Here we model the expansion of the QGP using viscous 
hydrodynamics with the bulk viscosity <^ set to zero. It is 
known that C scales as c^ — 1/3, so from Fig. [T]it is clear 
that for the LSM C ^ (see also Refs. 28, 29]). Recent 
Lattice QCD results indicate that the quark-gluon matter 
EOS departures from the conformal description |70l - [74| . 
although it is common practice to neglect ( in the hydro- 
dynamic equations while using an EOS inspired on Lat- 
tice QCD. An estimate for the temperature-dependence 
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and the relative importance of ^ can be obtained from 
(e — 3p)/T^, which measures the deviation from the con- 
formal limit. Fig. [7] shows this quantity calculated in 
the LSM as a function of temperature, for g — 3, 3.2, 3.4. 
It is seen that, for the range of temperatures of interest 
T > 130 MeV, (e - 3p)/T^ is small except at T < T^ 
where it rises abruptely reaching significant values. The 
increase in (e — 3p)/r* near Tc is sharper for larger values 
of g. It has been shown before that the hydrodynamic 
evolution is affected by bulk viscosity, for example reduc- 
ing the elliptic flow, inducing the phenomenon of cavi- 
tation or modifying the freeze-out process [ij, |7^, [S^l • 
Although it is reasonable to assume that these changes 
will not affect in a dramatic way the differences in V2 
and W4/(w2)^ that result from including or not the chi- 
ral fields as sources for the evolution of the quark fiuid, 
a more thorough investigation of this particular issue is 
necessary. 

We now go over to the calculation of the shear viscosity. 
In order to calculate rj we adopt the results of Ref [23] , 
which are based on the linearized Boltzmann equation 
in the relaxation time approximation. In the relaxation 
time approximation the shear viscosity is given by 



'n = -^ 



At r (Pp 
ST./ (27r)3^ 



P 



/o(l - /o) 



(A5) 



where r = t{T) is the collision time. t{T) is calcu- 
lated from the averaged cross sections a for quark-quark 
and quark-antiquark scattering processes including 1/A^c 
next to leading order corrections as: 



' = 6/o( 



^uu-^uu \ ^ud—>ud I ^uu- 



ud^-ud 



(A6) 



We refer the reader to Refs. [28|, |9J, |95| for details on 
the calculation of the a' s - we note that the chiral model 
used in these references is not very different from the 
LSM, the exchanges of a and tt mesons are modeled by 
contact terms. Note also that the cross sections are tem- 
perature dependent, not only because of phase space, 
but also because they depend on the constituent quark 
masses, whose temperature dependence is given by the 
LSM of the present paper. See also Ref. [231 fo'^ ^ related 
approach applied to the calculation of -q/s in the LSM of 
a pion gas. 



Appendix B: Dependence of results on cut— off for 

■q/s 

It was noted before that any hydrodynamical descrip- 
tion of matter is bound to become inapplicable at large 
viscosity. Since our simulations are purely hydrodynamic 
and the calculated rj/ s increases rapidly with decreas- 
ing temperature, a cut-off for the value of rj/ s must be 
imposed. The comparison of the results obtained from 
hydrodynamic and kinetic simulations to data provide a 



guide as to what is the value at which one should im- 
pose the cut-off. In our simulations we set the cut-off 
at -q/s = 0.4, which is in the range of values found to 
correspond to the breakdown of viscous hydrodynam- 
ics [50, 173). In this Appendix we show the results for 
hadronic observables calculated using different values of 
this cut-off, namely rj/ s < 0.3, 0.4 and 0.5, and analyse 
the dependence of the results on this choice for the cut- 
off. For simplicity, the results shown in this appendix 
correspond to g = 3.2. 
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FIG. 8: (Color online) V2 for charged-hadrons calculated tak- 
ing or not taking into account the chiral fields as sources for 
the hydrodynamic equations calculated with different choices 
for the cut-off value rj/ s = 0.3, 0.4, 0.5. The value of the 
coupling constant is 5 = 3.2. 
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FIG. 9: (Color online) «4/(«2)^ for charged-hadrons calcu- 
lated taking or not taking into account the chiral fields as 
sources for the hydrodynamic equations calculated with dif- 
ferent choices for the cut-off value rj/ s = 0.3, 0.4, 0.5. The 
value of the coupling constant is <; = 3.2. 

Figs. [5] and IH] show the charged-hadron V2 and 
''^4/(f^2)^, respectively, as a function of transverse momen- 
tum, obtained with cut-off values of rj/ s = 0.3, 0.4 and 
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0.5. It is seen that although the values of 172 and Vi/iv^f' 
change with the value of the cut-off, the relation between 
the results obtained including or not the chiral fields as 
sources in the hydrodynamic equations remains practi- 
cally the same. Consequently, the conclusions extracted 
from these results, which are discussed in the main text, 
do not depend on the precise value of the cut-off imposed 
on 7y/s, provided that g ^ 3.2 corresponding to a smooth 
crossover. 

As a final side remark, it is interesting to note that if 
one would attempt to extract the value of rj/s by match- 
ing V2 to data (with the remark mentioned in the In- 
troduction), the uncertainty in the extracted rj/s coming 
from the choice of cut-off in the range [0.3, 0.5] would be 
^ 50% with respect to the average value of rj/s. This 
relatively large value of the uncertainty in rj/s associated 



with the imposed cut-off clearly displays the sensitiv- 
ity of observables on the shear viscosity of the hadronic 
stage in collisions at y/ s^N = 200 GeV, a fact already 
discussed in Refs. |8l-[ll|. pointing to the conclusion that 
further study is needed to extract the precise temperature 
dependence of 77/s from heavy-ion collisions experiments. 
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